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X : Abstract 



In this paper, we investigate how a sizeable width-to-mass ratio for a gluino, 
as is for example realized in GMSB scenarios, could affect the discovery 
potential of gluinos at the LHC. More importantly, the influence of the 
gluino being "fat" on the standard mass and spin determination methods at 
the LHC are investigated. For this purpose, we focus on gluino production at 
the LHC, where we do not factorize the first step in the gluino decay cascade, 
but treat the following decay cascades step in factorization, including full 
spin correlations. The effects of sizeable width-to-mass ratios from a few up 
to 15-20 per cent on the endpoint of several mass determination methods as 
well as on means for discrimination between BSM spin paradigms like SUSY 
and UED are studied. 



^ juergen . reuterOdesy . de 
^daniel . wieslerOdesy . de 



1 Introduction and Finite Width Effects 



The Standard Model (SM) of particle physics describes all data in the field since today with a very 
good accuracy. Nevertheless, there are many reasons for physics beyond the SM (BSM), namely the 
non-existence of a dark matter particle within the SM, the insufficent amount of CP violation to explain 
the baryon-antibaryon asymmetry in the universe, and the instability of mass terms of fundamental 
scalar particles against radiative corrections. Supersymmetry has been one of the most favorable 
candidates to cure these problems. One of its main predictions to be tested best at a hadron collider 
is the existence of strongly interacting supersymmetric partner particles of the quarks and the gluon, 
namely squarks and the gluino. Up to now, in the runs at 7 and 8 TeV center-of-mass energy, the 
Large Hadron Collider (LHC) has not found any traces of supersymmetric particles, which means that 
these particles cannot be too light. Supersymmetry is part of a large class of BSM models, where new 
particles around the TeV scale have only weak interactions, like in models with extra dimensions, or 
Little Higgs models. These kind of models seems to be favored by electroweak precision data, as there 
are no (big) deviations or inconsistencies in the fit of the electroweak data to the SM. Generically, new 
particles in weakly interacting models show up as narrow resonances, where in most cases their width 
is below the detector resolution of the LHC experiments. This is different from strongly interacting 
scenarios like technicolor, composite Higgs models or conformal sectors where new particles or more 
like broad resonances (e.g. similar to the p resonance) or even similar to continuum-like excitations 
like in QCD or condensed matter physics. 

In the minimal supersymmetric SM (MSSM) and its most simple extensions, all new particles are 
consequently rather narrow resonances with a width-to-mass ratio, 7 := T/M, of the order of half a 
per cent or less. Besides the heavy Higgs bosons in certain parts of parameter space, the gluino g is 
the only particle that can get a sizeable width-to-mass ratio. In the case, the gluino is heavier than 
the squarks, there are many decay channels open for the gluino, such that the gluino can easily access 
a width-to-mass ratio of several, and even up to 15-20 per cent. The theoretical upper limit for this 
ratio (taking the gluino mass to inifinity or equivalently assuming quasi-massless squarks) is 32 per 
cent. More details under which conditions gluinos can become "fat" are summarized below, when we 
discuss the model setup for this study. 

Since the most severe bottleneck of simulations with multi-particle final states (at least, but not 
only at tree level) is the integration over the high-dimensional phase space, almost all SUSY studies 
have been performed with a factorized approach. This is motivated by the fact that a production 
process of two SUSY particles is followed by subsequent two- or three-body decays. The easiest 
aprroximation relies on the narrow-width approximation (NWA) which precisely does this factoriza- 
tion. This can be approved by folding in momentum smearing according to a Breit-Wigner propagator 
N{q)/{q'^ — m?' + imV), where q is the four momentum, m the mass and F the width of the associated 
intermediate particle. In many cases, the most severe effects come from the numerator factor N{q) 
which depends upon the quantum numbers of the intermediate state, and essentially contains the spin 
density matrix connecting the production and decay matrix elements. Spin and polarisation informa- 
tion for the intermediate particle are important e.g. in top or W decays, or in BSM physics for the 
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discrimination between different assignments of particle spins, e.g. between SUSY and UED. 

When the decay matrix element (squared) is integrated out, the small width T of the intermediate 
state particle is neglected, as well as the (usually) very small dependence of production and decay 
matrix elements on the momentum of the intermediate particle q. In a previous study [T] it was shown if 
one of the five following cases is met, the NWA cannot be applied: (1) the obvious violation, when 7 = 
T/M is not small, e.g. for the p resonance in QCD or a corresponding resonance in composite models, 
(2) if the masses of mother and daughter particle are near-degenerate, m ~ M (decay threshold), 
e.g. typically in UED, (3) if the mass of the intermediate is close to the kinematical limit of the 
experiment, M ~ ^/s (particle threshold), (4) if there are big interferences between different exclusive 
decay channels (quasi-combinatorial background), and finally (5) if there is a propagator non-separable 
from the matrix element (e.g. in the case of a non-trivial spectral density of particles or unparticles). 
In the SM, there are only a few examples where deviations from the NWA really matter, as many 
(fundamental) particles are relatively narrow (the W with 7 =2.5 per cent has the largest width-to- 
mass ratio) in the SM. However, even here sometimes the NWA had to be given up in order to get to 
predictions that are precise enough to match experimental analyses, e.g. e~^e~ — t- W^W~ — )• 4/ [2]. 
In general BSM scenarios, such cases appear much more often, cf. e.g. |3|4j . Though there have been 
systematic improvements of the NWA |5|6] up to the order 0{T/M), it is always preferrable to use 
full matrix elements for the processes under investigation. 

In this paper, we study whether and how much finite-width effects of a heavy gluino do affect its 
detectability, and its mass and spin determination at the LHC. Sec. [2] defines our benchmark SUSY 
models, the studied processes and the setup of our simulation. The results are shown in Sec. [3] for the 
effects of finite widths on the mass determination and in Sec. H] for the spin determination. Finally, 
we conclude in Sec. O 

2 Simulation Setup and Benchmark Model 

For the analyses of off-shell effects on endpoint and shape measurements of invariant mass distributions, 
we assume a SUSY scenario where there are considerable branching ratios of the gluino both into a 
first- or second-generation squark accompanied by a jet. Furthermore, we assume the presence of a 
so-called golden chain with subsequent decays of the squark into the second-lightest neutralino, which 
further decays into slepton and lepton, where the slepton then ends up in another lepton and the 
lightest neutralino. This demands for the specific mass hierarchy in the SUSY particle spectrum, 
niq^ > > '"^[^ > "^xi' which however is not uncommon in phenomenological SUSY models 

studied in the literature (early cMSSM scenarios such as SPSla [7] inspired the presence of this type 
of hierarchy). For the study of gluino width effects, we overlay the two main production processes 
of gluino pair production with the much more abundant - in the case of light (er) squarks - one of 
associated gluino-squark production: 




pp am + X 



(1) 
(2) 
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In the first case, Eq. ([I])), we take asymmetric decay chains, where one of the two signal gluinos decays 
into two hght (without loss of generality down) quarks and the lightest neutralino, while the other 
one decays into two bottom quarks and the second-to-lightest neutralino, which further decays via an 
intermediate (right handed) slepton to two corresponding leptons and a lightest neutralino: 

gi bbi bhxl bbl^ll bbl^Pxi (3) 
92 ^ ddL ^ ddxl (4) 

Notice the index i at the bottom squark owing to the fact, that we include both decay modes of the 
gluino into bi and 62- In the second case of squark-gluino associated production (Eq. ([2])), we simulate 
the prompt squark decay into a quark and the lightest neutralino 

qi/R qxi- (5) 

With the focus on this particular exclusive final state being the same for both production processes we 
omit additional complications from combinatorial ambiguities. It allows us to study the consequences 
of off-shell effects without having to suffer from SUSY or combinatorial backgrounds (these have 
been studied recently in the context of purely hadronic decay modes in |8|11| . Preliminary results on 
distortions due to off-shell effects have been shown in [9], while distortions of parton-level distributions 
due to exotic non-standard SUSY particles which could mimic off-shell and combinatorial effects have 
been studied in [lO] . 

The last decay steps in the cascade of equation, Eq. ([3|), are particularly well-known: they are 
identical to the golden chain with the replacement of a (first or second generation) squark by a 
sbottom: — ?• bi. This exclusive final state allows for a quite generic study of most of the mass 
determination methods in the literature |12j , while in principle simultaneously reducing combinatorial 
mis-assignments due to the possibility of 6-tagging. Furthermore, in the analysis of off-shell gluino 
width effects we are able to use methods based on the existence of sbottoms in the cascade. 

In the following, we introduce a benchmark scenario that allows for the abovementioned signal 
decay chain while capturing most of the relevant phenomenological features like a Higgs mass around 
125 GeV, heavy colored states beyond the LHC limits [13] and a rich phenomenology of decay patterns. 
To be as generic as possible, we decided to make use of the phenomenological MSSM with 19 free 
parameters (plDMSSM), without any high-scale relations among the parameters. All parameters 
determining the model are given at the electroweak/TeV scale. Their explicit values for the parameter 
point chosen here are given in Table [TJ The mass hierarchy for the decay chain above is valid for 
this parameter point, and the branching ratios of the four successive two body decay steps are of 
considerable size: 



g bbi 10 % 
g 662 07 % 



bi bxl 16 % 



62 ^bxl 34 



/o 



Xl ^e^el 42% 
Xl ^f^^fil 42% 



e| ^ e±xi 100 % 
^/^^X? 100% 



Given these values, the total branching fraction of our exclusive final state from signal gluinos decaying 
through the benchmark cascade above is roughly 7 %. This number has to be taken with a grain of 
salt, as the underlying concept of factorized cross sections for exclusive final states into cross sections 
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Table 1: Model parameters of the pMSSM parameter point studied in this paper. All numbers (except 
tan (3 ) are in units of Ge V. 

times branching ratio as within the NWA might no longer be a good approximation in the presence 
of off-shell contributions. However, in this study we investigate the effects arising from precisely 
these kind of contributions far away from the resonant pole of the propagator. A complete treatment 
should (at least) take into account fully differential four- (or more) particle final states. Since we use 
these figure merely as rough estimates to determine the actual number of events we expect from our 
exclusive decay cascade final state, we refrain from such a calculation. All masses of the spectrum were 
calculated using SOFTSUSY [T3], while the particle decay widths were obtained with SUSYHIT [15]. An 
overview of all model parameters is given in Table [TJ Feeding the obtained SLHA file |16|17j for the 
given mass spectrum into PROSPINO [18], we calculated the relevant cross sections for LHC at 14 TeV: 
for squark-gluino associated production 376.0 fb, for gluino pair production 47.6 fb, and neutralino- 
gluino associated production 3.7 fb, respectively. The sum we multiply by the total branching fraction 
of our exclusive final state (7 %) to arrive at an event number of roughly 9,000 for an overall integrated 
luminosity of 300 fb~^. Several diluting effects such as detector acceptance, b-tagging efficiencies and 
event selection criteria may further reduce this figure. For that reason, we took a conservative estimate 
and chose to analyse an event number of 5,000. The following simulations are based on the most recent 
version of the event generator WHIZARD [19]. All generated events were exported to HepMC [20] format 
and afterwards passed through a C++ analysis framework based on ROOT |21] . which was specifically 
developed and tailored to the investigation of deviations from off-shell contributions in spin and mass 
determination methods. 

As we want to investigate effects that come directly from off-shell propagators and interference 
effects, we compare full matrix-element calculations with factorized approaches in the NWA at parton 
level, in order to be able to disentangle their effects from pollution that stems from QCD radiation, 
hadronization and detector effects. To study a "fat" gluino, we use full matrix elements for the 
production including the first decay of the gluino. More precisely, the first part of Eq. ([3]) is completely 
calculated in one step, 

pp ^ (66x2) + (g/q) , (6) 

including all interferences. On the other hand, the successive decays of the second-to-lightest neu- 
tralino and the "spectator" gluino or squark are factorised with full spin correlations using the NWA. 
Furthermore, the phenomenological width-to-mass ratio 7 = F/M is scanned over the following values 

7 G {0.5%, 2.5%, 5.0%, 10.0%, 15.0%, 20%, 25%}. (7) 
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This set is chosen so as to resemble a broad range of possible widths that might be realised in nature. 
While the lowest value is a typical width for standard parameter points, values above 10 % appear 
when the gluino is clearly heavier than the squarks as in GMSB setups. There is also a particular 
phenomenological reasoning behind the range of values, which will be discussed at the beginning of 
the next section. Using this relative width ratio 7 we study the impact of off-shell contributions on a 
selected choice of mass and spin determination methods listed in the next two sections. 



3 Effects on Mass Measurements 

We begin our investigation with the study of several mass measurement variables. A large width F in 
the gluino propagator affects the momenta of both the intermediate sbottom and the near (bottom) 
quark 6„, as both are directly resulting out of the gluino propagator. The far bottom quark bf, 
coming from the sbottom decay, on the other hand is expected to receive only a minor modification 
and thus should not be part of distorted invariant mass distributions. Fig. [1] depicts the transverse 
momentum distributions of both the near and the far bottom quark and their distortion with respect 
due to different values of 7. The black (solid), red (short-dashed), green (dotted), blue (short-dashed- 
dotted), yellow (long-dashed-dotted), magenta (long-dashed-double-dotted) and cyan (long-dashed) 
line correspond to 7 = 0.5 %, 2.5 %, 5.0 %, 10.0 %, 15.0 %, 20.0 % and 25.0 %, respectively. 
As expected, the far bottom quark, is almost unaffected by a fat gluino as it originates from the 
subsequent decay of the sbottom. On the other hand, we observe that the near bottom quark exhibits 
a considerable distortion in its pT distributions, already visible by eye. While there is an obvious 
tendency for a severely increasing distortion within the first five values of 7 (up to 15 %), the two 
largest effective widths (20 and 25 %) exhibit only a slight further increase. These values are anyhow to 
some extent academic, since in realistic scenarios they are quite hard to realize: either these scenarios 
are already excluded due to the Higgs boson discovery at 125 GeV, or the gluinos are so heavy that 
they are out of reach of a 14 TeV LHC. In the first plot for the transverse momentum distributions, 
we included them to illustrate the effect of a moderate saturation. This observation motivates us to 
leave out the two highest values 7 = 20.0%, 25.0 % and from now on investigate the reduced range of 
values up to 15 %, also in order not to make the plots too crowded. 

As the far quark/jet is (almost) not affected by the off-shell effects, we restrict our investigation 
to that subset of mass determination observables introduced in [12] that contain the near quark bn- 
{ruhb, mb„i^iou,,rnb^e,high^ rnbbe,iow, rubbiMgh, rubbee}- These are the invariant mass distributions of 

either two b jets, one b jet and a lepton (with the softer or the harder jet, respectively), of two bs with 
a lepton, or finally two jets and two leptons. All these variables are tailor-made for exclusive decay 
cascades, which are shown in the next subsection, while the subsection after that shows the effect of 
broad gluinos on inclusive variables (like e.g. M^j). 



3.1 Exclusive Cascades 



rribb 



The first mass edge we study is the classical dijet endpoint of the first gluino decay step. 
Due to the scalar propagator of the intermediate bottom squark, we expect to see no effects of spin 
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Figure 1: Transverse momentum distributions of the two bottom quarks in the signal cascasde (near 
from the gluino and far from the shottom decay). The black (solid), red (short- dashed), green (dotted), 
blue (short- dashed- dotted), yellow (long- dashed- dotted), magenta (long- dashed- double- dotted) and cyan 
(long-dashed) line correspond to = 0.5 %, 2.5 %, 5.0 %, 10.0 %, 15.0 %, 20.0 % and 25.0 %, 
respectively. 

correlation, such that the shape of the distribution should resemble the well-known triangular nature 
of the di-lepton edge with a linear rise from to m^^, where a sharp cutoff marks the endpoint. 
Figure [2] shows on the left hand side the simulated distributions for five different values of 7, where 

"ibb Endpoint distribution 



800 p 




Figure 2: Left hand side: Invariant dijet mass of two bottom quarks, mj,i,, from the gluino decay 
for different values of 7 (line styles identical to Fig. Qp. On the right hand side, the distribution of 
endpoints obtained as fit parameters of the edge-to-bump method. 

the line colors are identical to Fig. [TJ One notices that the solid black line deviates only very little 
from the ideal triangular shape dictated by phase space (one basically has a small smearing at the edge 
of the distribution), which is a direct consequence of the relatively small width: F = 0.005 M. The 
situation drastically changes when the width is incrementally increased: the distortion steadily grows 
and starts to wash out the clear-cut edge structure with increasing 7. At high values of 7 (> 5 %) the 
distribution acquires an irreducible tail, which mimics other distorting effects such as combinatorial 
mis-assignments or detector smearing already at this early parton level stage. However, we want to 
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stress the fact, that this kind of deviation from off-shell contributions is irreducible in the sense, that 
it cannot be tuned away e.g. by other methods designed to minimize combinatorial problems. 

In the next step, it is our goal to quantify these intrinsic contributions by estimating the distortion 
of the shape in the vicinity of the endpoint as a function of the effective width-to-mass ratio 7. For 
this study, we use the edge-to-bump method [22] . It translates the human bias in fitting different lines 
to kinematical distribution and hence a systematic uncertainty into a statistical uncertainty over a 
sample of line/edge fit, thereby scanning of a variety of different fit ranges. It allows us to extract the 
edge information and to discriminate different endpoint behaviours in an unbiased way. In a nutshell, 
the approach fits a naive linear kink function O (1000) times and returns bumps at the most likely 
positions of inks (supposedly physical edges) in the original distribution. 

As a first estimate of the impact of width effects on the measurement, we consider the shift of the 
actual endpoint position with respect to the value of 7. Table [2] discloses in the second column these 
values obtained with our own implementation of the method described in detail in [9]. The according 
distributions are in Fig. [2] on the right hand side. While for small width-to-mass ratios of 0.5% the 
obtained value (708 GeV) is close to the template one (679.6 GeV), the endpoint positions for large 
of 7 are off by more than 200 GeV (920 GeV), using the same method and settings. 

An alternative, but related measure for the endpoint smearing is given by the size of the corre- 
sponding error estimates. Their increase with respect to 7 reflects the observation that the spread of 
endpoint values in the right plot of Fig. [2] is considerably enhanced for an increased effective width. 
While for 7 = 0.5% the purely statistical error is small, the sheer growth of the standard deviation 
going up to 7 = 15.0% by nearly two orders of magnitude serves as another good indication for a huge 
endpoint smearing. This raises an important point: the overly high confidence expressed through the 
small errors of endpoints for low values of 7 is a mere binning effect and does not represent a real- 
istic error estimate for sophisticated endpoint measurements. Moreover, these error estimates of the 
edge-to-bump method are purely statistical and refiect the transformation of a statistical uncertainty 
on an endpoint position onto the particular position of a mean value in a distribution of fit results. 
Since the binning of a histogram is bounded by experimental resolutions, a relatively large minimum 
bin size results in a systematic underestimation of the given errors. 

The usage of the endpoint position as a measure for distortion is based on just one fit parameter 
(the so-called p4, for details cf. |22|9| ) of the edge-to-bump method. A detailed observation of the 
edge fitting function suggests to make use of the sampling point parameters p2 and ps, which are the 
two linear slopes left and right of a possible edge or kink. Notwithstanding the fact, that the off-shell 
contributions tend to wash out the sharp edge and lead to smoother and longer tails, we propose to use 
the absolute difference of the two slopes as well as the ratio of the two slope parameters as measures 
to quantify the amount of edge distortion: 

Sd-=\P3-P2\ (8) 
Sr ■■= \P2/P3\ (9) 

For each edge fit, both values are calculated from the parameters returned by the edge-to-bump 
method. Ideally, for a pure phase space distribution of triangular shape and no smearing beyond the 
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sharp cutoff, the first slope is infinite and the second slope zero, hence the slope difference maximal 
(namely infinite) and the slope ratio minimal (namely zero). In that sense, the difference in slopes 
measures the sharpness of a kink in the distribution whereas the ratio returns information about the 
size of the second slope relative to the first one (pronouncedness of the kink). A ratio close to zero 
may thus be attributed to a tail-less distribution such as the triangular shaped one. Keep in mind, 
that due to the sanity checks of our method, \p2\ < \p3\ and hence Sr G (0, 1). The full treatment of 
large width effects on the other hand introduces a tail and smears the endpoint behavior, which gives 
rise to considerably smaller slope differences and higher slope ratios. If the returned fit value of such 
a slope difference is compatible with zero, the underlying distribution apparently lacks robust kinky 
features. Extraction of such shallow endpoints is therefore a very delicate task. 

In the last two columns of Tab. [2] we collect results for the slope differences and slope ratios of the 
invariant di-bottom mass for all five scan values of 7. In Fig. [3] we show the impact of the different 
values of 7: while we observe slope differences well above one and small slope ratios compatible with 
zero for small 7, at an effective width of already 5 % the mean Sd is reduced to a value smaller than 
one. The slope ratio exhibits an equivalent behavior: the mean value is increased by a factor of 10. 



7[%] 


^bb 


Sd 




0.5 


708.5 ± 0.9 


5.70 ± 2.90 


0.014 ± 0.008 


2.5 


740.9 ± 2.5 


1.71 ± 1.17 


0.064 ± 0.028 


5.0 


835.7 ± 19.2 


0.78 ± 0.24 


0.084 ± 0.023 


10.0 


886.5 ± 13.0 


0.67 ± 0.08 


0.141 ± 0.045 


15.0 


921.3 ± 25.9 


0.62 ± 0.04 


0.131 ± 0.035 



Table 2: Adapted mean values of endpoint positions (in GeV), slope differences (in 1/GeV) and slope 
ratios for invariant di-bottom mass mi,i,. 
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Figure 3: Statistically distributed slope differences (left) and ratios (right) for the dijet invariant mass 
ruhb obtained as fit parameters with the edge-to-bump method. 
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mb„e I Next, we investigate the mass edges obtained by the minimization and maximization over two 
possible lepton combinations with the near bottom quark. Their definitions are: 



rrib^eMgh = max [mb„^+, m^^^-] 



(10) 
(11) 



while their endpoints given are given by 
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By the very nature of the intermediate particles, there is a small correlation of the leptons and the 
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Figure 4: Left: Invariant mass of near bottom quark and lepton as a function of nibi^iow (upper plot) 
and rribi^high (lower plot) for different values of^. The black (solid), red (short- dashed), green (dotted), 
blue (short- dashed- dotted), and yellow (long- dashed- dotted) correspond to j = 0.5 %, 2.5 %, 5.0 %, 
10.0 %, and 15.0 %, respectively. Right: Distribution of the corresponding endpoints obtained as fit 
parameters with the edge-to-bump method for rrihi^iow (upper plot) and nihi^high (lower plot). 

bottom quark, that manifests itself in the shape of the distributions. Hence, with no pure phase 
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space shape, we expect the distortion to be less pronounced in comparison to the invariant di-bottom 
mass. Especially mh^j^io^ should not contain too much excess events in the upper parts, since these 
contributions will mostly be omitted due to the minimization procedure. The variable rrih^i^highj on 
the other hand, will severely be affected for the very same reason. The distributions in Fig. |4]on the 
left hand side (upper and lower plot) confirm these conjectures. The extracted endpoints depicted on 
the corresponding right hand sides of the upper and lower line in Fig. [4] and their numerical values in 
Table [3] also support this statement: while for ?7i;,„£^;ou) the overall endpoint variation with respect to 
7 is about 40 GeV, ^n-b^^^hig/i suffers from more than twice the endpoint shift with a value of about 80 
GeV. Comparing this to the expected ideal edge positions 



^bnl,l0W 



364.4 GeV 
493.1 GeV 



(14) 
(15) 



calculated from Eq. (jlip . we find that the discrepancy for m^^^ is indeed small, and the determined 
value is mostly in agreement with the expectation. For mh^i^high, however, the endpoint shift is 
as large as 15%. The adapted slope differences Sd show a similar picture: while the numbers of 
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Figure 5: Statistically distributed slope differences (left) and slope ratios (right) for rrihi^iow (upper 
plot) and rrihi^high (lower plot), obtained as fit parameters by means of the edge-to-bump method. 

the minimized distributions exhibit only a mild overall decline, the maximized version illustrates an 
apparent tendency of Sd to decrease with rising 7. As for the slope ratios Sr, the same holds, but to 
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an even greater extent, in that for nrihi^iow the three intermediate ratios are all in the same ballpark, 
whereas for mi,£^high, Sr increases by nearly one order of magnitude. 



7[%] 


"''b£,low 


Sd 




'"'be,hiah 






0.5 


341.4 ± 6.7 


3.25 ± 0.38 


0.014 ± 0.007 


499.8 ± 0.3 


2.35 ± 0.03 


0.017 ± 0.005 


2.5 


356.6 ± 2.7 


2.58 ± 0.47 


0.055 ± 0.025 


511.2 ± 3.6 


1.80 ± 0.26 


0.055 ± 0.018 


5.0 


355.5 ± 1.8 


2.82 ± 0.25 


0.041 ± 0.011 


527.0 ± 11.9 


1.55 ± 0.25 


0.054 ± 0.017 


10.0 


383.0 ± 39.6 


1.59 ± 0.77 


0.066 ± 0.024 


556.3 ± 13.4 


1.00 ± 0.44 


0.103 ± 0.048 


15.0 


377.2 ± 22.3 


1.47 ± 0.74 


0.124 ± 0.069 


579.0 ± 47.9 


0.79 ± 0.40 


0.158 ± 0.059 



Table 3: Adapted mean values of endpoint positions (in GeV), slope differences (in 1/GeV) and slope 
ratios for rriM^iow and rriM^high, respectively. 



rribbi I Extending the low- and high-type invariant masses from above by the additional (far) bottom 
quark, we arrive at a three-particle invariant mass, that has a similar feature as the ones just discussed: 
an intermediate neutralino propagator communicating spin correlations and allowing for a similarly 
altered shape compared to the triangular phase space. If we now start to gradually increase the 
effective width parameter 7, the off-shell contributions start to enter the game in a more severe way 
as for the two-particle invariant masses, rribe^iow and mhi^high- This can be understood in terms of 
the inclusion of nearest neighbors: since the invariant di-bottom mass is heavily distorted by width 
effects, and it is always included in both mbbe,iow and ruhbi^high, we in fact expect to observe large 
deviations. In that sense - and in contrast to mbi^iow ~ the minimization procedure over two possible 
lepton combinations is not able to suppress the appearance of these intrinsic contributions, as can 
be seen in Fig. [6l Regarding the shift of edge positions, the situation is comparable to the case of 
rribb'- both mbbi^iow and rubbi^high exhibit displacements of up to 150 and 180 GeV, respectively. For 
the highest value of 7 = 0.15 in Table HI we notice an exceptionally small error estimate for the 
endpoint position of mbbi^iow in contrast to an unusually large one for mbbi^high- The reasoning behind 
this is a trivial matter of statistics as is evident from Figure [6l while the lower of the two invariant 
masses has a sharp drop at the bin corresponding to the endpoint position, the higher distribution has 
two such fake kinks at around 1100 and 1250 GeV, respectively. These are purely statistical issues 
happening by chance and attributed to the low overall number of events of 5,000. Hence the gross 
under- and overestimation of the error estimates for both edges. The mean value however still captures 
the important feature of endpoint translation: a shift of the returned mean value of the edge-to-bump 
method of up to 180 GeV. Comparing these with the theoretically expected values 

= 868.6 (16) 
m'S>f%^gh = 996.6 (17) 

we again find agreement for small widths, but large deviations for high values of 7. In Fig. [7] and 
Table m the slope differences Sd (ratios Sr) of the two variables both display a clear decreasing (increas- 
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'bbl,low 



Endpoint distribution 




Figure 6: Left: Invariant mass of near bottom quark and leptons, varied as defined by mhM,iow (top) 
and mfjfji^high (bottom) for different values of ^ (line style identical to Fig. Right: Distribution of 
the corresponding endpoints obtained as fit parameters with the edge-to-bump method for mbu,iow (top) 
and mbu^high (bottom). 



7[%] 


'"'bbl,low 


Sd 










0.5 


845.2 ± 11.6 


2.87 ± 0.50 


0.025 ± 0.011 


1009.7 ± 00.4 


3.47 ± 0.10 


0.009 ± 0.005 


2.5 


864.4 ± 02.2 


1.79 ± 0.91 


0.046 ± 0.020 


1017.0 ± 01.7 


2.11 ± 0.77 


0.041 ± 0.012 


5.0 


900.8 ± 23.5 


1.45 ± 0.22 


0.067 ± 0.022 


1060.0 ± 10.6 


1.72 ± 0.32 


0.050 ± 0.012 


10.0 


979.4 ± 34.8 


0.62 ± 0.29 


0.156 ± 0.051 


1120.8 ± 11.3 


0.94 ± 0.22 


0.126 ± 0.031 


15.0 


1002.3 ± 4.9 


0.71 ± 0.07 


0.163 ± 0.041 


1184.3 ± 62.9 


0.63 ± 0.30 


0.175 ± 0.052 



Table 4: Adapted mean values of endpoint positions (in GeV), slope differences (in 1/GeV) and slope 
ratios for mm,iow and rribbiMah- 

ing) trend for growing (falling) effective width factors 7 and thus confirm our choice as a parameter 
quantifying the size of the off-shell effects on the quality of the mass determination. 



The last three-particle invariant mass we investigate is the combination of the near bottom 
quark and the lepton pair. Without the particular need for distinction between the two leptons, we 
expect to observe only a moderate distortion for large gluino widths, which might be weakened by 
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Figure 7: Statistically distributed slope differences (left) and slope ratios (right) for mhbi^iow (top) and 
f^hu^high (bottom), obtained as fit parameters of the edge-to-hump method. 

the fact that the first hnear slope of the undistorted distribution is not maximally steep since mb^u is 
not built out of direct next neighbors in the decay chain. An example of a maximally steep endpoint 
behavior is the undistorted invariant di-bottom mass, whose first slope parameter would ideally be 
infinite due to the triangular (phase-space) shape. Fig. [8] depict in the upper line mi,^u and the 
corresponding endpoint distributions, numerical values of which are given on the left-hand side of 
Table [H Despite a fairly accurate endpoint estimate for up to 7 = 2.5% (in comparison with the 
theoretical value of 578.8 GeV), there is still a deviation of substantial size for large widths of up to 
100 GeV. The two slope parameters both behave similarly to what we have already seen with other 
invariant mass variables (cf. Fig. [9]): an average slope difference Sd of well above two is reduced to a 
value below one, that is compatible with zero within two standard deviations. For small widths, the 
slope ratio Sr is also close to zero, but increasing 7 results in a steady growth of s,. of more than one 
order of magnitude in size to values of up to 16 %. 



Combining all objects from one cascade, we obtain the variable m^bee, which has a very clear- 
cut endpoint and a steep edge structure, since it is constructed out of next neighbors only. The 
sharper the drop of the distribution at the vicinity of the edge, the larger we expect the impact of 
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Figure 8: Left: Invariant mass of bottom quark(s) and leptons: mh^u (top) and nibbu (bottom) for dif- 
ferent values of J. The black (solid), red (short-dashed), green (dotted), blue (short- dashed- dotted), and 
yellow (long- dashed- dotted) correspond to ^ = 0.5 %, 2.5 %, 5.0 %, 10.0 %, and 15.0 %, respectively. 
Right: Distribution of the corresponding endpoints obtained as fit parameters of the edge-to-bump 
method for mb„u (top) and rubbu (bottom). 

off-shell contributions to be, even so if just one out of four propagators, in our case the signal gluino, 
is affected. From the lower line of Fig. [8l the amount of distortion is already visible by eye. A more 
quantitative statement is given in terms of numerical values on the right-hand side of Table [5j While 
the estimated endpoints for small values of 7 are in gross agreement with the theoretically expected 
value of 1092.7 GeV, the largest deviation of m^^^ for 7 = 15% is up to 172 GeV. A similar picture 
as for mb„u is found for the slope parameters (cf. also Fig. [9|): ranging from just below five down to 
well below one, the slope difference exhibits an even larger spread of values. The slope ratios on one 
hand are compatible with zero for the smallest width 7 = 0.5%, but, on the other hand, increase by 
a factor of 15 for the largest off-shell contribution. 

3.2 Inclusive Approaches 

Up to now, we only studied exclusive invariant mass variables, i.e. combinations of objects arising 
from one cascade side from a decay of just one single mother. In our case, this was taken to be the 
gluino, which has a large variety of possible decay patterns and consequently a plethora of interesting 
invariant mass combinations. Restricting the analysis to just one particular decay cascade entails 
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7 [%] 


^bll 






^bbll 






U.o 


5d5.4 ± 1.7 


o o'v 1 r» 1 c! 

Z.J7 ± U.15 


0.010 ± O.OOz 


1079.3 ± 2.4 


A C A 1 n 1 

4.64 ± 0.51 


A A1 1 1 r\ Ann 

0.011 ± 0.009 


2.5 


568.5 ± 3.4 


1.92 ± 0.08 


0.080 ± 0.037 


1105.1 ± 1.3 


2.22 ± 1.22 


0.035 ± 0.013 


5.0 


594.5 ± 2.0 


1.56 ± 0.04 


0.050 ± 0.012 


1150.2 ± 19.9 


1.34 ± 0.70 


0.078 ± 0.020 


10.0 


640.0 ± 12.0 


0.94 ± 0.25 


0.142 ± 0.064 


1193.7 ± 37.4 


1.00 ± 0.44 


0.163 ± 0.041 


15.0 


668.6 ± 40.5 


0.65 ± 0.33 


0.157 ± 0.050 


1251.2 ± 30.9 


0.78 ± 0.17 


0.182 ± 0.036 



Table 5: Adapted mean values of endpoint positions (in GeV), slope differences (in 1/GeV) and slope 
ratios for rrihu and rubbu ■ 
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Figure 9: Statistically distributed slope differences (left) and ratios (right) for rribu (top) and m^hu 
(bottom) obtained as fit parameters of the edge-to-bump method. 

several problems, the largest is presumably combinatorics. Consider for example the symmetric case 
of two identical decay chains of the type we have analyzed so far, where not one, but two gluinos decay 
into two bottom quarks, two leptons and the lightest neutralino. All variables we have just discussed 
assume that a differentiation between the two cascades is somehow given (e.g. the partitioning of four 
leptons into 2x2 leptons). However, a priori, there is no general recipe that always allows for such 
a correct assignment. This severely affects the usability of these exclusive variables and is known as 
the combinatorial problem. 
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Figure 10: Mt2 distribution (upper left) and its statistically distributed endpoints (upper right), slope 
differences (lower left) and slope ratios (lower right) obtained as fit parameters of the edge-to-bump 
method for five different values of ^. The color scheme is the same as in Fig. 0. 



An alternative, less exclusive approach is Mj'2 as introduced in [23], which is a grouping of 
the measured objects into two subsets and a following minimazation procedure over the invariant 
transverse masses. Although M'j'2 is similarly affected by the combinatorial problem (since here the 
visible momentum has to be split into two separate sides as well) there are methods which address 
this issujl]. In the remainder of this chapter we circumvent these kind of combinatorics through two 
non-identical decay chains, which allow us to directly concentrate on the off-shell effects. Hence, the 
partition 
given by 



partition of visible momenta into two sides for the application of Mr2(p|,2)Plis)|'T' "^x) trivially 



p^2 = {b,b,i^,n (18) 
P^2 = {q,q} , (19) 

where the light-flavor partonic jets are assigned to one side and the remaining objects to the other. 
The input test mass of the invisible sparticle was set to = 150 GeV close to the true value 
of the benchmark scenario. Figure [10] illustrates in the upper left plot the impact of the effective 
width parameter 7 on the distribution of Mt2- similar to the exclusive invariant mass variables 



^e.g. MtGoh [21], which is the minimum of A/t2 for all possible momentum assignments into two partitions 
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discussed above, relative widths larger than 2.5 % result in substantial smearing and a long tail in the 
distributions, which is absent at lower values. Owing to the general minimization procedure of Mt2, 
the effect is somewhat smaller than for exclusive variables such as mi^f, or m^fe^^, which is underlined 
by the numerical results given in Table [H The endpoint M^2^ is given by the parent sparticle mass, 



7[%] 




Sd 




0.5 


1293.7 ± 0.8 


5.92 ± 0.21 


0.009 ± 0.003 


2.5 


1315.0 ± 10.8 


4.42 ± 0.54 


0.029 ± 0.017 


5.0 


1381.0 ± 27.4 


1.82 ± 1.03 


0.054 ± 0.027 


10.0 


1526.0 ± 18.6 


0.70 ± 0.20 


0.072 ± 0.029 


15.0 


1589.5 ± 31.4 


0.87 ± 0.40 


0.079 ± 0.022 



Table 6: Adapted mean values of endpoint positions (in GeV), slope differences (in 1/GeV) and slope 
ratios for Mt2 ■ 

thus we expect to observe a clear edge structure at the gluino mass of 1277 GeV. Up to a slight 
overshoot, this is in gross agreement with the values given for 7 < 5%. However, effective widths of 
up to 15 % lead to an edge shift of nearly 300 GeV. This is the largest translation of an endpoint 
we have obtained so far and can be understood in terms of the underlying topology: in this inclusive 
scenario both partitions are affected by off-shell contributions through non-resonant parts of the gluino 
propagators, and hence they have their share in the distribution by means of the definition of Mj'2- 
To this extent, it was necessary to create a special event sub-sample, which slightly differs from the 
one introduced in the beginning. With respect to Eq. ([6]), the factorization of the spectator gluino 
was exchanged with by full matrix-element calculation according to 

PP ^ (bbx^) + (qQX?) , (20) 

since the light quarks now also play a role in the construction of Mt2- The six-particle final state 
(with subsequent factorized decay of the second-to-lightest neutralino X2) necessary to be able to 
fully analyze all off-shell effects in this inclusive scenario. The slope parameters in Table [6] confirm 
the strength of the deviation: while the differences exhibit a significant drop for rising values of 7, 
the ratios Sr show a minor but steady increase. 

In a first summary, mass determination methods generally all suffer from off-shell effects. The 
amount of smearing in the distribution depends on the exclusiveness or inclusiveness of the variable. In 
general, the more particles originating from one or even two fat gluinos are involved in the construction 
of the mass variable, the more distortion or smearing appears in the distributions. The errors in the 
mass determination already at parton-level can be as large as 10-15 per cent. 

4 Effects on Spin Determinations 

After a possible discovery of any new physics beyond the SM, the next steps in determining the under- 
lying model characteristics are the measurements of masses and spins of novel particles. While many 
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mass determinations rely on endpoint positions of invariant mass distributions, the nature of the un- 
derlying spin is encoded in the shape of those (and especially angular) distributions and as such more 
delicate to differentiate. Hence, after discussing the effects of off-shell contributions on mass determi- 
nation variables, in this section we turn to spin measurements and carefully analyze the effects of a fat 
gluino onto several methods designed to distinguish a hypothetical SUSY signal from an equivalent 
one of Universal Extra Dimensions (UED). In general, spin studies compare shapes of distributions 
by choosing one particular type of mass spectrum which is either of UED or supersymmetric nature. 
Though typical UED spectra are far more compressed resulting in softer decay products and thus 
require for more comprehensive analyses (cf. e.g. |25|26] ). we choose to stick to hierarchical SUSY like 
spectra to get a comprehensive comparison. Moreover, as we want to emphasize the difference of spin 
in intermediate propagators and their impact on invariant mass distributions before comparing it to 
contributions arising from off-shell effects, we construct a particular UED model, which inherits all 
masses and width parameters from our SUSY benchmark model and hence allows us to use the partic- 
ularly interesting decay chain already known from the mass measurement section with the following 
replacements: {g,bi,X2,lR,X2) — ^ {g^^\b^^\ Z^^\l^^\'j^^^). Thus, in this cascade, edges of invariant 
masses stay the same but shapes thereof are expected to drastically change. In contrast to the gluino, 
the Kaluza-Klein (KK) gluon will retain a small effective width of 7 = 0.5 % throughout this analysis. 
As for the technical side, we use an adapted version of the minimal UED model [27], implemented into 
WHIZARD using the FEYNRULES [28] interface for WHIZARD [29]. In the sequel, we start to analyze shape 
asymmetries based on exclusive invariant quark-lepton masses before turning to hadronic correlations 
inside a single cascade and finally investigating the impact of non-resonant contributions on inclusive 
angular distributions. 

4.1 Shape Asymmetries 

Many studies of spin measurements rely on the specific decay topology of the golden chain and make 
inherent use of invariant mass shapes therein as discussed above. Since we are interested in effects 
emerging from a gluino, we concentrate on an extended version: 

g gnQL Qnqjxl Inqj^^it qnqft^t^xi ■ (21) 

On the basis of this decay chain, we analyze to what extent the spin determination methods proposed 
in |26] are affected by the off-shell contributions from a fat gluino. The approaches studied within that 
paper were designed to discriminate signatures of a supersymmetric gluino from the ones emerging 
from a Kaluza-Klein gluon excitation in models of UED. 



A^(mf,£) At first, we investigate the special bottom-lepton asymmetry with the definition: 



^±^^^^^^<Mpil^±Ip^ . (22) 
da/d'mf^l+ + da/dmi,^~ 

We follow the assumptions in [26], namely that the bottom (instead of anti-bottom) quarks are uniquely 
identified through a lepton charge tag of the 6-tagging algorithm. Due to the Majorana nature of the 
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Figure 11: Invariant masses m^i+ (upper left) and mf,£- (upper right) with a bottom quark and their 
asymmetry for 5k (lower left) and 25k (lower right) events. The black (solid) and blue (dotted) 
line correspond to ^ = 0.5 % and 15.0 %, respectively. The UED sample is given by the red (dashed) 
line. 

gluino, all bottom quarks are near bottom quarks in 50 % of all decays, i.e. they are produced in the 
first two-body decay step. Hence, visible effects from off-shell contributions of a fat gluino propagator 
are expected to influence half of the invariant mass shapes. Fig. [TT] depicts both the invariant mass 
distributions of ?7i;,^+ and m^(^- as well as the bottom-lepton asymmetries, for 5k and 25k events. To 
keep the plots clear and at the same time display all relevant information, we refrain from using all 
widths, but rather restrict ourselves to the most extreme values of 7 = 0.5 and 15%, given by the 
black (solid) and blue (dotted) lines as well as the UED sample, depicted by the red (dashed) line. 

To clarify whether the deviation of the large-width sample is merely a statistical effect, we arti- 
ficially increased the event number by a factor of 5 up to 25k. As it turns out in Fig. [11] (on the 
lower right) there is indeed a subtle effect that is observable in the large width SUSY sample, which 
is well beyond the size of fluctuations, although both of the two invariant mass distributions show no 
strong discrepancies with respect to the two different values of 7. Nonetheless, the minimum plateau 
from 300 to 400 GeV as well as the subsequent rise from 500 to 800 GeV are both reduced by up to 
one third in magnitude. As is evident through direct comparison with the superimposed UED sample 
however, it is obvious that large off-shell contributions are not endangering a possible discrimination 
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Figure 12: Invariant masses mfj^£+ (upper left) and m^^^- (upper right) with the softer of the two 
bottom quarks and their asymmetry for 5k (lower left) and 25k (lower right) events. The color 
scheme corresponds to the one in Fig. [71] 



of the fundamentally different spin scenarios, using this specific variable. 



A^{mi,i) Next, we investigate the impact another bottom-lepton asymmetry. The difference to the 
first asymmetry discussed above is given by the spectrum dependent property that the softer b-quark 
may coincide with the nearer b-quark. The definition of this asymmetry is given by: 

A±(^ ^ _ da/dmf,^i+ - da/dm,,^^- 

K^bsl) — , / , , , , , ■ (2c5J 

da / dmi,^i+ + da / drufj^i- 

In fact, in our scenario, this assignment is true most of the time: for small width 7 = 0.5% in 
4557/5000 ~ 90% of all events, the near bottom quark is also the softer one. Apparently, this changes 
when the width is increased, as is illustrated in Fig.[TJ However, even for the largest width of 7 = 15%, 
in 4109/5000 ~ 80% of all events this assumption is still correct. Obviously, this characteristic will 
drastically change when the mass gap between the gluino and the sbottom is increased. The nearer 
quark will become harder due to a larger phase space in the gluino decay, and the value of will 
consequently be reduced. In the upper line of Fig. [12] we show both invariant masses, mj,^£+ and 
mi,^£-, which exhibit the typical smearing behavior for large widths encountered more often in the 
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mass measurement section above. Already by eye, the asymmetries in the lower line of Fig. [12] have 
a comparable pattern: both invariant mass distributions have similar shapes and the corresponding 
asymmetries are thus rather small. Regardless of the size of the effective width 7, the distortion of 
mi,_^£+ mimics the one of mf^^i-, and the same holds for the asymmetries. Consequently, the small 
event sample of 5k is fully compatible with a vanishing asymmetry throughout the complete range of 
the histogram, not only for the two SUSY samples, but also for UED. Moreover, the larger samples 
of 25k events also have only minor deviations from = 0, which might be attributed to statistical 
fluctuations. Altogether, we find the size of this second bottom-lepton asymmetry to be of negligible 
size compared to the already small deviations found in [26]. Hence, it is not surprising that a steady 
distortion for both mi,^i+ and mf,^^- results in a negligible change of an already very small asymmetry. 
We conclude that this last asymmetry is neither preferable in terms of discriminative power between 
SUSY and UED nor for suffering from smearing due to off-shell effects. After all, this should be 
attributed to the specific kind of underlying mass spectrum. 

4.2 Hadronic Angular Correlations 

It was further proposed to analyse purely hadronic correlations such as the average pseudo-rapidity 
fjhi) or the difference of azimuthal angles A;/);,;, of the two bottom quarks [26]. The first of these 




Figure 13: Average pseudo-rapidity (left) and azimuthal distance (right) of both bottom quarks for 5k 
(top) and 25k (bottom) events. The coloring corresponds to the one in Fig. [771 
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observables exhibits no visible distortion with respect to the maximal width 7 in the supersymmetric 
sample, although a small difference to the UED case is perceptible. We confirm the results of [26], 
where a slightly more central behavior of the two 6-quarks was found for SUSY signals. For the 
azimuthal distance, we find a non-negligible deviation (cf. Fig. [T3l) of the large width SUSY sample 
from the standard SUSY sample. In addition, the former turns out to resemble the shape of UED 
events. As a next step, for each of the two variables the following asymmetries are defined |26j : 

N{\f],,\ < 1) - iV(|%fc| > 1) 



N{\fibb\ < 1) + N{\fibb\ > 1) 



(24) 
(25) 



N{A^bb < + N{A<Pbb > 7T/2) ' 
which were proposed to obtain an additional measure that allows for a discrimination between a 
standard SUSY and a UED signal. We apply these variables to our three scenarios, namely the 
standard and off-shell SUSY as well as the UED samples. Numerical values for all of these are 
given in Table [7] for 5k and 25k events, respectively. The value of the asymmetry of the average 



sample 


5k 25k 


A^ (std) 
At (ofs) 
At (ued) 


0.627 ± 0.017 0.628 ± 0.008 
0.645 ± 0.017 0.645 ± 0.008 
0.567 ± 0.016 0.557 ± 0.007 


A^ (std) 
At (ofs) 
At (ued) 


0.014 ± 0.014 0.005 ± 0.006 
-0.047 ± 0.014 -0.052 ± 0.006 
-0.042 ± 0.014 -0.039 ± 0.006 


At (std) 
At (ofs) 
At (ued) 


0.194 ± 0.015 0.180 ± 0.007 
0.125 ± 0.014 0.129 ± 0.006 
0.003 ± 0.014 0.008 ± 0.006 



Table 7: Numerical figures for the three asymmetries At, At and At defined in equations 
and (j28p for different scenarios: std, ofs and ued correspond to the standard = 0.5%J and off-shell 
(7 = 15% ) SUSY as well as UED event samples. Errors are purely statistical. 

pseudo-rapidity fjhi, only exhibits a marginal increase for the off-shell SUSY sample compared to the 
undistorted standard one. Both coincide within at most two sigma of the purely statistical error. The 
UED case on the other hand has less central values than the off-shell SUSY parts. This difference to 
both the standard SUSY sample as well as to the UED sample should be attributed to the fact that the 
off-shell contributions tend to harden (at least the first) decay product(s) (recall the pT distributions 
of the near and far bottom quarks in Figure [1]), and hence allow for even more central values of r]. 

The situation looks a lot less promising for the azimuthal distance A(f)bb of the two bottom quarks. 
Here, we find that off-shell contributions in the large-width SUSY sample drive the asymmetry At 
such that a discrimination between SUSY and UED is no longer possible. Although the shapes of the 
distributions differ only by a moderate amount visible by eye in all three scenarios (cf. Fig. [T3l upper 
right), the off-shell SUSY sample is washed out in such a way so as to drive the quantitative numerical 
estimate of the asymmetry negative by the same amount as in the case of UED (cf. Table [7|). 
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Figure 14: cos6* for 5k (left) and 50k (right) events. The color scheme is the same as in Fig. 



4.3 Inclusive Angular Distributions 

Finally, we investigate angular correlations of the initially produced particles, which in our case corre- 
sponds to gluinos or KK gluons. Although we are not able to reconstruct the complete mother particle 
momenta due to missing energy, the first emitted partons of each decay cascade should still possess 
an observable angular correlation among each other. This was first used in the variable 



cos Oil = tanh 



(26) 



in a study of slepton pair production [30], and later adapted to general colored SUSY production [3T] . 
The adapted method was applied to fully hadronized inclusive signal event samples with gluino and 
squark contributions or the corresponding equivalent for UED, where the largest discriminative power 
is attributed to the squark and KK quark signatures. Our study on the other hand is based on a 
parton- level analysis, and we restrict ourselves to subsamples with gluinos, since we aim to assess the 
impact of their off-shell contributions on the method. More precisely, we apply the variable to the 
exclusive gluino benchmark process in our scenario introduced in the beginning in the following way: 



cos e*gg 



tanh 



(27) 



where qi = min(6i,62) and qr = min(gi,g2) are the softer of the two quarks from each cascade side. 
While in |31j the largest contributions arise from prompt squark decays to quarks and the lightest 
neutralinos, and thus cos 9*^ is chosen to be applied to the two hardest objects, we make particular use 
of the (unfortunately spectrum-dependent) approach of selecting the softer quarks to be attributed 
to the first (near) gluino or KK-gluon decay products. These are furthermore assumed to inherit 
features of the initially produced mother particles. Fig. [T3] illustrates the behavior of cos 9*^ for the 
two SUSY scenarios with different widths (black and blue) and the UED scenario (red) for 5k and 
50k events, respectively. We included the high statistics sample to show the asymptotic behavior of 
the fundamentally different models. The distortion due to off-shell contributions is apparent, and 
although they tend to wash out the distribution to less central values, the differentiation with respect 
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to UED is not endangered in the exclusive gluino subsample. A quantification of this statement can 
be obtained through the definition of the following asymmetry |31| : 

^ jV(|cosg;,| <0.5)-iV(|cosg*,| >0.5) 

Af(|cos^*J < 0.5) +iV(|cos0*J > 0.5) ' ^ 

whose values for the three cases are given in Table [71 The initial observations from Fig. [T3] are 
confirmed: while the value of is indeed reduced by one third in the standard SUSY sample with 
respect to the off-shell SUSY sample, the UED case is compatible with a value of zero. Hence we 
conclude, that although an there is an apparent modification of cos 9*^ due to finite- width effects, it 
is not threatening the discrimination of SUSY and UED models. 



5 Conclusions 

In this paper, the impact of off-shellor finite- width effects from a broad ( "fat" ) gluino on a represen- 
tative selection of observables for mass and spin determination is investigated with a discrete set of 
effective width-to-mass ratios 7 = F/M E {0.5%, 2.5%, 5.0%, 10.0%, 15.0%}. Such larger widths of 
gluinos are not uncommon in certain regions of parameter space, and are e.g. quite natural in GMSB 
scenarios. The basis for this analysis are SUSY signals from gluinos where production and the first 
decay step are simulated with a full matrix-element calculation to account for non-resonant contri- 
butions in the gluino propagator and interferences, distorting kinematic observables. All subsequent 
decays are then factorized using the NWA, including full spin correlations. An adapted edge-to-bump 
method is employed to quantify deformations in kinematic distributions, which arise in scenarios with 
large effective width-to-mass ratios 7. In this approach, several parameters from a naive linear kink 
function fitted to the corresponding variables are extracted and utilized to quantify errors in the mass 
and spin determination arising from finite-width effects. In general, mass measurement observables are 
much more obviously affected in that the endpoint smearing dominates the distributions of otherwise 
sharp edge structures already at moderate width-to-mass ratios of 7 = 5%, and steadily increases 
throughout higher values of 7. As a result, the endpoint positions obtained with the edge-to-bump 
method are drastically shifted already at the parton level. The difference and the ratio of the two 
slope parameters in the fit are additional measures of distortion which similarly indicate a washed-out 
endpoint behavior. 

Methods of spin determination are on the other hand less affected by propagator contributions far 
off the mass-shell. Shape asymmetries as well as angular correlations exhibit only small deviations 
of event samples with broad gluinos from the ones with a narrow resonance. In the special case of 
azimuthal distance, the difference of these two event samples however is of considerable size, and 
although the shapes of the broad SUSY as well as the UED sample still differ, the numerical values 
of the associated asymmetry coincide within their purely statistical errors. 

In summary, our studies show that broader BSM resonances as in our case a fat gluino lead to 
severe distortions of kinematic distributions which are the basis for many mass and spin determination 
methods. These deviations from non-resonant contributions arise at the fundamental parton level in 
simulations where production and decay are simulated with the full matrix element. Many existing 
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studies on the other hand utihse the NWA, whose predictions vigorously differ from the full calculation. 
Hence a correct treatment of such effects is of crucial importance for both mass and spin measurements 
so as to not misinterpret the potential signals of much sought-after new physics. While the spin 
determination might not suffer too much in such cases, broader particles like a fat gluino might make 
life hard for mass determinations at the LHC. 
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